Classical and Quantum Chaos and Control of Heat Flow 



O 
O 

(N 

» 

o 
O 
o 



Giulio Casati 1,2,3 and Carlos Meji'a-Monasterio 4 
1 Center for Nonlinear and Complex Systems, Universita degli Studi dell'Insubria, Como, Italy 
2 CNR-INFM and Istituto Nazionale di Fisica Nucleare, Sezione di Milano 
i Department of Physics, National University of Singapore, Singapore 117542, Republic of Singapore and 
4 Dipartimento di Matematica, Politecnico di Torino, 
Corso Duca degli Abruzzi 24 1-10129 Torino, Italy 
(Dated: February 6, 2008) 

We discuss the problem of heat conduction in classical and quantum low dimensional systems 
from a microscopic point of view. At the classical level we provide convincing numerical evidence 
for the validity of Fourier law of heat conduction in linear mixing systems, i. e. in systems without 
exponential instability. At the quantum level, where motion is characterized by the lack of expo- 
nential dynamical instability, we show that the validity of Fourier law is in direct relation with the 
onset of quantum chaos. We then study the phenomenon of thermal rectification and briefly discuss 
the different types of microscopic mechanisms that lead to the rectification of heat flow. The control 
of heat conduction by nonlinearity opens the possibility to propose new devices such as a thermal 
rectifier. 

PACS numbers: 05.70.Ln,05.45.Mt,44.10.+i 
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I. INTRODUCTION 

The origin of the macroscopic phcnomenological 
"laws" of thermodynamic transport is still one of the 
major challenges to theoretical physics. In particular the 
issue of energy (heat) transport, in spite of having a long 
history, is not completely settled 0, @ • Given a partic- 
ular classical, many-body Hamiltonian system, neither 
phenomenological nor fundamental transport theory can 
predict whether or not this specific Hamiltonian system 
yields an energy transport governed by the Fourier law 
J = — kVT, relating the macroscopic heat flux to the 
temperature gradient VT Heat flow is universally 
presumed to obey a simple diffusion equation which can 
be regarded as the continuum limit of a discrete random 
walk. In consequence, transport theory requires that the 
underlying deterministic dynamics yield a truly random 
process. Therefore, it is not mere idle curiosity to won- 
der what class, if any, of many-body systems satisfy the 
necessary stringent requirements. 

In spite of intense investi gat ions in rece nt y ear s, [H, la , 
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the precise conditions that a dynamical system of inter- 
acting particles in ID must satisfy in order to obey the 
Fourier law of heat conduction are still not known. How- 
ever, the general picture which emerges is that, for sys- 
tems with no globally conserved quantities (i.e., globally 
ergodic), positive Lyapunov exponent is a sufficient con- 
dition to ensure Fourier heat law. This appears quite 
natural. 

Indeed modern ergodic theory tells us that for K- 
systems, a sequence of measurements with finite preci- 
sion mimics a truly random sequence and therefore these 
systems appear precisely those deterministically random 
systems tacitly required by transport theory. In fact 
the thermal conductivity has been studied for a Lorentz 
channel -a quasi one dimensional billiard with circular 



scatterers- and it was shown to obey Fourier law |15j . 
Yet we do not have rigorous results and in spite of several 
efforts, the connection between Lyapunov exponents, cor- 
relations decay, diffusive and transport properties is still 
not completely clear. For example, in a recent paper |23j 
a model has been presented which has zero Lyapunov ex- 
ponent and yet it exhibits unbounded Gaussian diffusive 
behavior. Since diffusive behavior is at the root of nor- 
mal heat transport then the above result [2^ constitutes 
a strong suggestion that normal heat conduction can take 
place even without the strong requirement of exponential 
instability. Moreover, the models in 0, |2^| are nonin- 
teracting and thus, the condition of Local Thermal Equi- 
librium (LTE) is not satisfied. Noninteracting particle 
systems are certainly less realistic. They are good mod- 
els for gases in the Knudsen limit in which the mean free 
path is larger than the characteristic length of the con- 
tainer and their use in other physical situations is open 
to criticism [24j. Therefore, it is interesting to investi- 
gate up to what extend one can simplify the microscopic 
dynamics, and yet obtain a normal transport behavior. 

An additional interesting problem, almost completely 
unexplored, is the derivation of Fourier law from quan- 
tum dynamics. Indeed at the quantum level, inves- 
tigations have been mainly focused on linear response 
theory 25] . The possibility to derive the Fourier law by 
directly establishing the dependence of J on VT in a 
nonequilibrium steady state calls directly in question the 
issue of quantum chaos. In this connection we recall that 
a main feature of quantum motion is the lack of expo- 
nential dynamical instability 26], a property which is at 
the heart of classical dynamical chaos. 

The aim of this paper is to discuss the transport of heat 
in several classical and quantum microscopic models by 
direct numerical simulation of energy flow in the system 
in contact with thermal baths. At the classical level our 
results in section [H] provide convincing evidence that lin- 
ear mixing systems (i.e. with zero Lyapunov exponent) 
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obey Fourier law. At the quantum level we study in sec- 
tion [TO] an Ising chain of L coupled spins 1/2 in a tilted 
magnetic field. We provide clear numerical evidence of 
the validity of Fourier law at the onset of quantum chaos. 
In section lTVl we then investigate the possibility to control 
the energy transport and discuss different microscopic 
mechanical models in which thermal rectification can be 
observed. The possibility to control heat conduction by 
nonlinearity opens the way to design a thermal rectifier, 
i.e. a lattice that carries heat preferentially in one di- 
rection while it behaves as an insulator in the opposite 
direction. 



II. DYNAMICAL INSTABILITY AND FOURIER 

LAW 




FIG. 1: The geometry of the model. Particles move in the 
region outside the triangular scatterers. The x coordinate 
goes along the channel and y is perpendicular to it. The 
two heat reservoirs at temperatures Tl and Tr are indicated. 
The length of each cell is I = 3, the base of the triangles is 
a — 2.19 and the distance between the two parallel lines is 
d — 1.8. The geometry is then uniquely specified by assigning 
the internal angles 9 and (j> 



A Lorentz gas consists of noninteracting point particles 
that collide elastically with a set of circular scatterers on 
the plane. Motivated by the ergodicity and mixing prop- 
erties of the Lorentz gas, in [15| a channel geometry of 
this model was considered to study the problem of heat 
conduction. By imposing and external thermal gradient 
it was found that in a Lorentz channel the Fourier law 
is satisfied. Yet this system is not described by LTE. A 
modification of this model in which particles and scat- 
terers can exchange energy through their collisions ap- 
peared in 20] . This effective interaction leads to the 
establishment of LTE. As a consequence, this model has 
proven to reproduce realistically macroscopic transport 
in many different situations pi) . In particular, the valid- 
ity of Fourier law was verified, indicating that the inter- 
action among the different particles is not fundamental 
for the observation of normal transport. 

To clarify the role of dynamical instability on the valid- 
ity of Fourier law, let us consider a two dimensional bil- 
liard model which consists of two parallel lines of length 
L at distance d and a series of triangular scatterers (Fig. 
1). In this geometry, no particle can move between the 
two reservoirs without suffering elastic collisions with the 
triangles. Therefore this model is analogous to the one 
studied in [15| with triangles instead of discs and the es- 
sential difference is that in the triangular model discussed 
here the dynamical instability is linear and therefore the 
Lyapunov exponent is zero. Strong numerical evidence 
has been recently given 18] that the motion inside a tri- 
angular billiard, with all angles irrational with tt is mix- 
ing, without any time scale, (sec also 19]). Moreover, an 
area preserving map, which was derived as an approxi- 
mation of the boundary map for the irrational triangle, 
when considered on the cylinder shows a nice Gaussian 
diffusive behavior even though the Lyapunov exponent of 
the map is zero |23| ■ It is therefore reasonable to expect 
that the motion inside the irrational polygonal area of 
Fig. 1 is diffusive thus leading to normal conductivity. 

We consider a channel of total length L = Nl where N 
and I are the number and the length of the fundamental 
cells. For the irrational angles we take 6 = 1)tt/2 



and <f) = 1. By increasing L, the number of particles per 
cell must be kept constant. However, since the particles 
do not interact one may consider the motion of a single 
particle over long times and then rescale the flux. 

Heat conductivity of this model has been numerically 
studied in |28[. Heat baths have been simulated with 
stochastic kernels of Gaussian type, namely, the proba- 
bility distribution of velocities for particles coming out 
from the baths is 

™-¥"(4)-™-isr -*(-#) 

(1) 

for v x and v v , respectively. 

Since the energy changes only at collisions with the 
heat baths, the heat flux is given by 

1 Nc 

j(t c ) = -J2( AE h> ( 2 ) 

c fc=i 

where (AE)k = Em — E out is the change of energy at 
the fc-th collision with the heat bath and N c is the total 
number of such collisions which occur during time t c . We 
have checked that for sufficiently long integration times 
(> 10 10 time unit) the system reach a stationary value. 

In Fig. 2, we plot the heat flux J as a function of the 
system size N. For the case of irrational angles, the best 
fit gives J oc A~ 7 , with 7 = 0.99 ± 0.01. The coefficient 
of thermal conductivity is therefore independent on A, 
which means that the Fourier law is obeyed. 

As a consistency check we have also performed an inde- 
pendent check of the above results via a Green-Kubo type 
formalism)^] from which consistent results with those 
presented above have been obtained. 

A completely different behavior is obtained when the 
angles 9 and </> are rational multiples of tt. The case with 
6 = tt/5 and cj> = 7r/3, is shown in Fig. (triangles), from 
which a divergent behavior of the coefficient of thermal 
conductivity k ~ jV - 22 is observed, indicating the ab- 
sence of Fourier law. 
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FIG. 2: Scaling behavior of the stationary heat flux J as a 
function of the system size for the case of irrational angles 
of Fig. 1 (circles) and for the rational case (triangles) . TV is 
the number of fundamental cells. The particle density was set 
to 1 particle per cell independently of TV. The best-squares 
fitting gives a slope —0.99 ± 0.01 for the irrational case and 
—0.78 ± 0.01 for the rational one. 



models is that in case (ii) the total momentum is con- 
served while in case (i) it is not. We recall that in several 
recent papers 0, |3(], it has been suggested that to- 
tal momentum conservation does not allow Fourier law 
and this may explain the lack of Fourier law for the one 
dimensional alternating mass model. 

To clarify the role of total momentum conservation, in 
[32| a model which is identical to the alternate mass hard- 
point gas but without total momentum conservation was 
considered. Numerical results clearly indicate that this 
model, contrary to the translationally invariant model, 
obeys the Fourier law. 

In perspective, these results demonstrate that diffu- 
sive energy transport and Fourier law can take place in 
marginally stable (non-chaotic) interacting many-particle 
systems. As a consequence exponential instability (Lya- 
punov chaos) is not necessary for the establishment of the 
Fourier law. Furthermore, they show that breaking the 
total momentum conservation is crucial for the validity 
of Fourier law while, somehow surprisingly, a less impor- 
tant role seems to be played by the degree of dynamical 
chaos. 



In conclusion, when all angles are irrational multiples 
of 7T the model shown in Fig. 1 exhibits Fourier law of 
heat conduction together with nice diffusive properties. 
However, when all angles are rational multiple of tt, the 
model shows abnormal diffusion and the heat conduction 
does not follow the Fourier law. 

One may argue that the model considered here is some- 
how artificial and far from realistic physical models. In- 
deed, noninteracting particle systems are certainly less 
realistic as in general, LTE is not established. However 
the problem under discussion is quite delicate and con- 
troversial and our main purpose is to understand which 
dynamical properties are necessary and sufficient to de- 
rive Fourier law. 

In this respect we mention two different billiards like 
models that are genuinely interacting many-particle mod- 
els and that share some of the properties of the model 
discussed here. The first one is the Id hard-point par- 
ticles with alternating masses that consists of a one- 
dimensional chain of elastically colliding free particles 
with alternate masses m and M ElGi. For this model 
all Lyapunov exponents are zero like the irrational tri- 
angle channel. In | 27| the problem of heat conduction in 
this model was considered. It was found that in the alter- 
nating mass model the heat current follows a power-law 
behavior J cx TV~" with a ~ 0.745. In contrast with the 
irrational triangle channel, the alternating mass model 
does not obey Fourier law. 

Therefore we have two models which are both mixing 
and without exponential instability: (i) the triangular 
billiard channel |28| (see also |2j|), which exhibits Fourier 
law and (ii) the alternate mass hard point gas mo del 
in which the coefficient of thermal conductivity diverges 
with the system size. The difference between the two 



III. FOURIER LAW AND THE ONSET OF 
QUANTUM CHAOS 

In section [n] we have shown that strong, exponen- 
tial unstable, classical chaos is not necessary (actually, 
strictly speaking, is not even sufficient |£j|) for nor- 
mal transport. In this connection, a main feature of 
quantum motion is the lack of exponential dynamical 
instability 18]. Thus it is interesting to inquire if, and un- 
der what conditions, Fourier law emerges from the laws 
of quantum mechanics. 

We consider an Ising chain of L spins 1/2 with cou- 
pling constant Q subject to a uniform magnetic field 
h = (h x ,0,h z ), with open boundaries. The Hamiltonian 
reads 

L-2 L-l 
n=0 n=Q 

where the operators a n — (a*, <r^, <r*) are the Pauli ma- 
trices for the n-th spin, n = 0, 1, ... L — 1. We set the 
coupling constant Q — 2. In this system, the only triv- 
ial symmetry is a reflection symmetry, a n — aL—i—n- 
Moreover the direction of the magnetic field affects the 
qualitative behavior of the system: If h z = 0, the Hamil- 
tonian (j3J corresponds to the Ising chain in a transversal 
magnetic field. In this case the system is integrable as 
© can be mapped into a model of free fermions through 
standard Wigner- Jordan transformations. When h z is 
increased from zero, the system is no longer integrable 
and when h z is of the same order of h x quantum chaos 
sets in leading to a very complex structure of quantum 
states as well as to fluctuations in the spectrum that are 
statistically described by Random Matrix Theory (RMT) 
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FIG. 3: Nearest neighbor level spacing distribution P(s) for 
the chaotic (left), intermediate (center) and integrable (right) 
chains. P(s) was obtained by diagonalizing the Hamiltonian 
J3J for a chain of length L = 12, averaging over even and odd 
parity sub-spectra. The curves correspond to Avd (dashed 
line) and to Pp (dotted line). 



|33| . The system becomes again (nearly) integrable when 
h z ^> h x . Therefore, by choosing the direction of the ex- 
ternal field we can explore different regimes of quantum 
dynamics. 

The onset of quantum chaos is commonly studied in 
terms of the nearest neighbor level spacing distribution 
P(s) that gives the probability density to find two ad- 
jacent levels at a distance s. For an integrable system 
the distribution P(s) has typically a Poisson distribution 
Pp(s) = exp (— s). In contrast, in the quantum chaos 
regime (for Hamiltonians obeying time-reversal invari- 
ance), P(s) is given by the Gaussian Orthogonal En- 
semble of random matrices (GOE). In this case, P(s) 
is well-approximated by the Wigner surmise Pwr>(s) = 
(ns/2) exp(— 7rs 2 /4), exhibiting the so-called "level repul- 
sion" . 

In Fig. |3 we show the results of our numerical simula- 
tions of system © for three different values of the mag- 
netic field: (i) chaotic case h — (3.375, 0, 2) at which the 
distribution P(s) agrees with GOE and thus corresponds 
to the regime of quantum chaos, (ii) integrable case 
h = (3.375,0,0), at which P(s) is close to the Poisson 
distribution, and (Hi) intermediate case h = (7.875,0,2) 
at which the distribution P(s) shows a combination of 
(weak) level repulsion and an exponential tail. 

In order to study energy transport we need to couple 
both ends of the chain of spins to thermal reservoirs at 
different temperatures. In [3~i| we have devised a sim- 
ple way to simulate this coupling, namely the state of 
the spin in contact with the bath is statistically deter- 
mined by a Boltzmann distribution with parameter T. 
Our model for the reservoirs is analogous to the stochas- 
tic thermal reservoirs defined by Eqs.n an d we thus call 
it a quantum stochastic reservoir. In what follows we use 
units in which Planck and Boltzmann constants are set 
to unity h = &b = 1 • 

The dynamics of the spins is obtained from the unitary 
evolution operator U(t) = exp(— iTit). Additionally the 
leftmost and rightmost spins of the chain are coupled to 
quantum stochastic reservoirs at temperatures P^ 1 and 
/3~ respectively. 

The action of the quantum reservoirs is as follows: at 



times t = nr with n integer the boundary spins couple 
to the quantum reservoirs. A local measurement of the 
boundary spins is performed and their new state is set to 
the state down (up) with probability /j (1 — /x) where the 
probability fJ-((3) depends on the canonical temperature 
of each of the thermal reservoirs as 



-fjjh _|_ e f3jh 



(4) 



The value of r controls the strength of the coupling 
to the bath. In 34] it was shown that averages over the 
ensemble of quantum trajectories or time averages over a 
quantum trajectory are sufficient to reconstruct a ther- 
mal density matrix operator that correctly describes the 
internal thermal state of the system in and out of equi- 
librium. Moreover, the time averaged expectation values 
of the local energy density can be used as a consistent 
canonical local temperature. In order to compute the 
energy profile we write the Hamiltonian @ in terms of 
local energy density operators H n : 



H n - -Q°~n°~n4 



h 

2 ' {<Jn 



'n+lj 



(5) 



The local Hamiltonian H n (defined for < n < L — 2), 
gives the energy density between the n-th and (n+ l)-th 
spins. In terms of eq. (JSJ the Hamiltonian of the system 
can be rewritten as 



H 



L-2 



'-(ay 



(6) 



In Fig. 0] we show the energy profile (H n ) for an out of 
equilibrium simulation of the chaotic chain. In all non- 
equilibrium simulations, the temperatures of the baths 
were set to T\ = 5 and T r = 50. After an appropriate 
scaling the profiles for different sizes collapse to the same 
curve. More interesting, in the bulk of the chain the 
energy profile is in very good approximation linear. In 
contrast, we show that in the case of the integrable (inset 
I) and intermediate (inset II) chains, no energy gradient 
is created. 

We now define the local current operators through the 
equation of continuity: dtH n — i\H,H n ) — — (J n+1 — J„), 
requiring that J n = [H n ,H n _x\. Using eqs. JSJ and © 
the local current operators are explicitly given by 

Jn = h x Q (o£_! - <+i) l<n<L-2. (7) 

In Fig. [3] we plot J/ AT as a function of the size L of 
the system for sizes up to L — 20. The mean current 
J is calculated as an average of (J n ) over time and the 
L — 8 central values of n. Since AL = L — 8 is an ef- 
fective size of the truncated system, the observed 1/AL 
dependence confirms that the transport is normal. On 
the other hand, in integrable and intermediate chains we 
have observed that the average heat current does not de- 
pend on the size J cx L°, clearly indicating a ballistic 
transport. 
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FIG. 4: Out of equilibrium energy profile {H n ) for the chaotic 
chain. The temperatures of the baths were set to 71 = 5 and 
T r = 50. Results for chains of size L = 8 (crosses), L = 14 
(open circles) and L — 20 (solid circles) are shown. The 
dashed line was obtained from a linear fit of the data for 
L — 20 for the L — 4 central spins. Insets (I) and (II) show 
the energy profile for the integrable and intermediate cases 
respectively, for L — 15. 



FIG. 6: Temperature profile of the chain with parameters 
Dl ~ Dr — 0.5, N = 128, a = 1, K = 0.3, M = 8, and T L = 
0.16, Tr — 0.15. The Morse potential constants are Do — 0.5 
- homogeneous system- (circles), Do = 1.2 (diamonds) and 
Do = 0.8 (+) 



IV. THERMAL RECTIFICATION 
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FIG. 5: Size dependence of the energy current in the chaotic 
chain with Ti = 5 and T r = 50. The dashed line corresponds 
to 1/AL scaling. In the inset, the size dependence of the 
energy current is shown for the integrable (circles) and the 
intermediate (squares). 



In conclusion we have shown that Fourier law of heat 
conduction can be derived from the pure quantum dy- 
namical evolution without any additional statistical as- 
sumption. Our results suggest that onset of quantum 
chaos is required for the validity of Fourier law. 



We now turn our attention to the possibility to control 
heat flow. To this end we start with the model introduced 
in 1331 and consider the Hamiltonian 



H 



71=1, JV 



Pn 

2m 



V n (y„ 



1 



Vn-l) 



(8) 



which describes a chain of N particles with harmonic 
coupling of constant K and a Morse on-site potential 
Vn{yn) — D n [eT anVri — l) 2 . This model was introduced 
for DNA chains where m is the reduced mass of a base 
pair, y n denotes the stretching from equilibrium position 
of the hydrogen bonds connecting the two bases of the 
n-th pair and p n is its momentum |3t| . In the context of 
the present study, model JSJ can simply be viewed as a 
generic system of anharmonic coupled oscillators, the on- 
site potential arising from interactions with other parts 
of the system, not included in the model. The Morse po- 
tential is simply an example of a highly anharmonic soft 
potential, which has a frequency that decreases drasti- 
cally when the amplitude of the motion increases. 

We consider the out-of-equilibrium properties of model 
||SJ| by numerically simulating the dynamics of the TV par- 
ticle chain, coupled at the two ends, with thermal baths 
at different temperatures. We thermalize at Tl, and Tr 
the first and the last L particles by using Nose-Hoover 
thermostats chains .331! or a Langevin description when 
we investigate cases very far from equilibrium. The baths 
temperatures Tl , Tr are never large enough to drive the 
system beyond the thermal denaturation temperature T c 
above which the mean value of y n diverges |35| . 

We compute the temperature profile inside the sys- 
tem, i.e. the local temperature at site n defined as 
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T n = m(y n 2 ), where ( ) stands for temporal average, 
and the local heat flux J„ = K(y' n (y n — y n -i))- The sim- 
ulations are performed long enough to allow the system 
to reach a non-equilibrium steady state where the local 
heat flux is constant along the chain. 

As a preliminary step we have considered the homoge- 
neous case in which D n = D, a n = a, n = 1, ..N. Here, 
as expected, we have detected a temperature gradient 
inside the chain, and we have verified that the system 
obeys the Fourier law of heat conduction [3^ . 

For the heterogeneous case we divide the chain be- 
tween the thermostats in three regions in which D n 
takes different values. In the left and right regions, 
D n = D L for ri = 1, ... (TV — M)/2 and D n = D R for 
n = (N + M)/2 + 1, . . . , N, respectively. For the M sites 
of the central region we set D n = D$. For the whole 
chain a n = a. 

In Fig. 38], we show the temperature profiles for 
three different values of Dq. Clearly a conductor- 
insulator transition occurs in the heterogeneous system 
as a function of the inhomogcneity measured by Do. 
We have obtained that the averaged heat flux J, which 
is uniform along the whole chain in the steady state, 
decreases by two orders of magnitude when Dq — D 
(D = D L = Dn) increases from to 0.7. This transition 
is explained in terms of a linearized version of (JSJ 3^ : for 
the nonlinear Hamiltonian Eq. [SJ effective phonon bands 
can be defined whose size and position depend on the 
microscopic properties of the interaction (in particular of 
coefficient D) and on the temperature. Using this depen- 
dence in a smart way one can design a thermal rectifier 
by choosing appropriate parameters Dl, Dn, Dq and a, 
depending on the range of temperatures in question. We 
intend for a thermal rectifier a device in which, for a fixed 
external temperature gradient, the magnitude of the heat 
current depends on the sign of the gradient. 

For instance, for the heterogeneous model one can de- 
sign a thermal rectifier by setting a strongly nonlinear 
region sandwiched between two weakly anharmonic left 
and right domains. In the presence of a thermal gradi- 
ent in the central part, the effective phonon frequencies 
evolve in space in a way that depends on the orientation 
of the gradient. This can provide either a good match- 
ing of the bands at the interfaces, leading to a thermal 
conduction across the system, or a complete mismatch 
leading to poor conduction [33. We exemplify this sit- 
uation in Fig. \7\ where the temperature profiles for the 
two orientations of the thermal gradient are shown. For 
the parameters values of Fig. the heat current changes 
by a factor of about 2 when the direction of the thermal 
gradient is reversed. 

Recently we have shown that using different micro- 
scopic models for the interaction like e.g. a Frenkel- 
Kontorova on-site potential, it is possible to strongly im- 
prove the efficiency of the rectifier 39] . Even though 
the underlying microscopic mechanism discussed in [39j 
is different from the one in ^33 it is based on the same 
general idea: the temperature dependence of the phonon 
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FIG. 7: Temperature profiles in a "thermal rectifier" for 
two opposite orientations of the thermal gradient. The dash- 
dotted lines show the borders of the different regions in the 
lattice. The thermostats have a size L = 8, the central re- 
gion a size M — 128, and the left and right regions contain 
64 particles. The coupling constant is K = 0.18, and the 
parameters of the Morse potential are Dl = 4.5, o?l = 0.5, 
Do = 0.7, Q = 1.4, D R = 2.8, a R = 0.5 in the left, cen- 
tral and right regions respectively. The temperatures of the 
thermostats are T = 0.1 and 0.45. When the high temper- 
ature is on the right side of the lattice, the average flux is 
J = 0.146 10 -3 while, when the thermal gradient is reversed 
the flux drops to J = 0.755 10 -4 . Note in this case the dis- 
continuities at the interfaces between region which attest of 
the bad energy transfer at these points. 

band. This dependence is due to the nonlinearity of the 
potential and therefore it should be possible to observe 
the rectifying effect in any nonlinear lattice. 

But thermal rectification is not exclusive to nonlin- 
ear lattices. Recently we also have shown that this phe- 
nomenon can also be observed in billiard systems of inter- 
acting particles . This model is based on the effective 
interaction of a gas of particles with the scatterers of a 
Lorentz channel |2(J, |21| . We have shown that the inter- 
acting character of the particles gives rise to dynamical 
memory effects that depend on the local thermodynam- 
ical fields, in particular on the local temperature gradi- 
ents. These memory effects have been described in detail 
in [12, El. We have used this dependence to set up a par- 
ticular geometry for which large thermal rectifications 
(the heat currents for both orientations of the thermal 
gradient differ by several orders of magnitude) are ob- 
tained [40j] . The possibility to obtain large rectification 
of the heat flow in billiard systems has raised a great 
interest because they are more easily realizable experi- 
mentally in the rapidly expanding field of nanophysics. 



V. CONCLUSIONS 

We have discussed the problem of heat conduction in 
classical and quantum low dimensional systems in rela- 
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tion to the dynamical properties of the system. At the 
classical level we have provided convincing numerical evi- 
dence for the validity of Fourier law of heat conduction in 
linear mixing systems, i.e. in systems without exponen- 
tial instability. As a consequence exponential instability 
(Lyapunov chaos) is not necessary for the establishment 
of the Fourier law. Moreover, we have shown that break- 
ing the total momentum conservation is crucial for the 
validity of Fourier law while, somehow surprisingly, a less 
important role seems to be played by the degree of dy- 
namical chaos. 

At the quantum level, where the motion is charac- 
terized by the lack of exponential dynamical instability, 
we have shown that Fourier law of heat conduction can 
be derived from the pure quantum dynamical evolution 
without any additional statistical assumption. Similarly 
to our observations in classical models, our results for a 
chain of interacting spins suggest that in quantum me- 
chanics, which is characterized by the lack of exponential 
dynamical instability, the onset of quantum chaos is re- 



quired for the validity of Fourier law. 

We have also discussed the phenomenon of thermal 
rectification in different classical models and discussed 
different types of microscopic classical mechanisms that 
lead to the rectification of heat flow. Even though all 
these mechanism are rooted in one way or another in 
the nonlinearity of the microscopic dynamics, we believe 
that, as for the case of Fourier law, thermal rectification 
is also possible at the quantum level. Some steps in this 
direction have already be taken li^L El ] 
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